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I. INTRODUCTION 



The Standard Model (SM) of particle physics provides a successful description of all 
experimental high-energy data to date. However, despite its success, many fundamental 
questions remain unanswered, such as the origin of electroweak symmetry breaking, the 
nature of neutrino masses, the large hierarchy between the electroweak and the Planck 
scales and the origins of dark matter and the cosmo logical constant. Attempts to address 
these questions have led to a wide range of new physics theories, most of them predicting 
new phenomena at the TeV scale. The Large Hadron Collider (LHC) at CERN is currently 
probing this scale and will hopefully make it possible to discover, constrain and/or exclude 
some of the proposed theories beyond the Standard Model (BSM). 

The Higgs boson as the last piece of the SM has been the searched for over several 
decades now. After the unsuccessful search for it at the LEP e~^e~ collider as well as the 
Tevatron (though there are some hints in a slight excess now in a broad mass range from 
120-140 Gev) several models that explain electroweak symmetry breaking (EWSB) as well 
as unitarization of tree level scattering amplitudes have been constructed that either work 
without a scalar state being present or sport invisible or undetectable decays of those states. 
One type of these models are the so-called HEIDI models [l|-[3| (an onomatopoeic version 
of the acronym for higher dimensions), based on an older model presented in |5|. The basic 
idea of this class of models is that the Higgs field mixes with a certain number of inert 
(singlet) fields, the number of which can be increased until they become a continuum. This 
admixture with such a set or continuum of fields can distort the properties of the physical 
Higgs particle(s), i.e. its production cross sections, its decay width and branching ratios. 
At the moment, there seems to be some evidence for a SM-like Higgs boson around 125 
GeV, which could have modified couplings to SM particles, as the diphoton rate seems to 
be rather enhanced, while the WW* decay seems to be suppressed with respect to the SM. 
Recently, the aforementioned models have been modified in order to precisely explain these 
deviations from the SM values 

An interesting feature of this type of model is that it may hide the source of electroweak 
symmetry breaking from discovery at the LHC. We know that, if the electroweak symmetry 
is broken explicitly in the Lagrangian (rather than through an asymmetric vacuum), the 
scattering of longitudinal electroweak vector bosons grows quadratically with energy and 
exceeds the perturbative unitarity bound at just over a TeV. As a result, we know that 
some new physics must appear below that scale in order to unitarize the theory. The 
common assumption is that we cannot fail to discover this new physics at the LHC. HEIDI 
type models, however, have areas of parameter space which are consistent with precision 
constraints, but may have such a broad signal as to be undetectable at the LHC with 
current methods [l|, Q . 

In this work we will be examining a HEIDI type model which includes an infinite but 
discrete set of additional scalar resonances. The model is realized by means of a spatial and 
compact extra dimension with the SM being confined to a four-dimensional brane within 
the five- dimensional spacetime. The only field which is allowed to propagate into the bulk 
of the extra dimension is a singlet scalar field coupled to the scalar sector of the SM through 
a trilinear coupling. After Kaluza-Klein (KK) reduction, we obtain an infinite number 
of new scalar fields mixing with the SM Higgs boson. Depending on the point in the 
parameter space, many KK states could be well below the TeV scale, and the emerging 
collider phenomenology of the scalar sector can thus be very different from the SM. In 
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order to assess the impact of the modified scalar sector on the Higgs searches, we study the 
so-called golden channel gluon fusion process. 

The outline of this paper is as follows: in Sec. [TT] we present our class of HEIDI models 
including their Lagrangian and derive the mass eigenvalues and eigenstates of the scalar 
sector together with its couplings to the SM particles. We discuss the parameter space of 
HEIDI models and examine different characteristic incarnations of the scalar sector depend- 
ing on the hierarchies in the parameter space. The simulation of such models is potentiall 
very demanding as the number of physical particles in the spectrum can be very large and 
depends on the point in parameter space, and our way of tackling these technical difficulties 
is discusse in Sec. Illli In Sec. ITV] our results for the golden channel of a Higgs boson decaying 
into ZZ* into four leptons at the LHC are shown, which is (if the Higgs mass is not too close 
to the LEP limit) the cleanest and most background-free channel. Finally, we conclude. 



II. THE MODEL 



A. Lagrangian and spectrum 

The compact HEIDI model discussed in this paper is a representative of the class of 
models suggested in Ref. The model is a renormalizable extension of the Standard Model 
formulated on Minkowski spacetime plus an additional spatial extra dimension compactified 
on a circle with radius R and denoted in the following by the coordinate y. In this model 
the Standard Model is confined to the y = brane, with the Higgs sector supplemented by 
an additional scalar gauge-singlet field Q{x,y) which is the only field allowed to propagate 
in the bulk of the extra dimension. The action of the scalar sector is chosen^ to be 




with $ being the Standard Model Higgs doublet (for later convenience we use a non-canonical 
normalization of $) and S<p the action describing the SM Higgs sector. From Eq. ([T]), it is 
evident that the only difference with respect to the SM is the additional sector Sq, which 
contains a trilinear coupling g between the Higgs field and the new scalar singlet Q. 

In order to work out the equivalent four-dimensional theory, we decompose the new scalar 
field in terms of its KK excitations^ 

n{x, y) = y~^ r- ^k{x) cos —y , with = ^/{l + Sko)7T . (2) 



We omit any self-interactions of $7 in the same spirit as Ref. 1] as they would spoil renormalizability. 
^ The sine modes in the KK decomposition decouple from the SM, so we omit them. 
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Inserting Eq. ([2]) into Eq. ([T]), we obtain the four- dimensional action 

Sn = I d'x f2 (^^d,nkd>^nk - ^mlnl + _ ^$t$^^^ ^ (3) 



k=0 

where the masses of the KK modes are given by 
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Like in the SM, spontaneous symmetry breaking occurs if /x^ > (see Appendix lA II 
for a detailed discussion of the potential and the vacuum expectation values). Rotating to 
unitarity gauge and shifting the scalar fields by their vacuum expectation values (vev), we 
get 

Hx)^nit. = f ^ + ] > ^k{x) =Wk + ujk{x) , (4) 



where the vevs are given by 

V = , ^ I with a = — coth (mhRir) | , = J^^ . (5) 

Defining the Higgs mass in the usual way, 

mh = V2\v , (6) 

the four-dimensional Lagrangian for the scalar fields, expanded around the correct vacuum, 
reads 



+ (^d^^kd^^k - ^mlul + ^^^ huk + ^ fi^Uk ] . (7) 



In the equation above, we have omitted constant terms as well as all the terms involving 
gauge bosons included in the covariant derivative of Eq. ([1]) (which are unchanged with 
respect to the SM). Note that Eq. (|5]) implies a condition on the quartic Higgs coupling A, 
the trilinear coupling g"^, the bulk mass rrih and the compactification scale R~^, 

X > a = — coth (miyRTr) . (8) 

As shown in Appendix lA 1\ violation of Eq. (|8]) leads to an unstable theory with a potential 
which is not bounded from below^. 

Owing to the trilinear coupling between $ and Q in the original action, Eq. (j7]) contains a 
mixing between h and the Uk which therefore do not correspond to mass eigenstates. While 



^ With the exception X = a and < 0, in which case the potential is bounded and has a single minimum 
at V — Wk — 0- 
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this mixing is resummed to obtain an effective Higgs propagator in Ref. [l|, we here take 
a different approach and diagonahze the mass matrix exphcitly. This calculation leads to 
the mass eigenvalues as the zeroes of a transcendental function (see Appendix lA 21 for 
details), 

f{\l) =ml + -^£= cot (RnJxl - - = . (9) 

It is easy to see from monotony arguments that the spectrum consists of exactly one mode 
A^ in each interval 



' + R2 < K < < + ^ with A: > 1 , (10) 



as well as an additional mode An below mi. 



B. SM couplings and parameter space 

The mass eigenstates (pk belonging to the mass eigenvalues A^ are linear combinations of 
the fields h and Uk, 

oo 

ct>k = eoh+Y.^t^i-i . (11) 
1=1 

Using the hermiticity of the mass matrix, we can invert this relation to derive the relations 
between h and the in terms of the mass eigenstates 

oo oo 

h = Y,eo<Pi , o;, = 5^^1+10,. (12) 

1=0 1=0 

From Eq. f|T2|) . it is clear that the SM fields couple to the physical fields just like they would 
couple to the SM Higgs boson, but with a factor for each 0/ leg at the vertex. As the 
couplings among the scalar fields are slightly more complicated, we have moved the detailed 
Feynman rules to Appendix IA3I 

According to Eq. ([T]), Eq. and Eq. (jH]), the model can be parameterized by the SM 
parameters rrih and v together with the bulk mass m^, the trilinear coupling g"^ and the 
compactification scale R~^, which is the parameterization that we will adopt for the rest of 
this work. However, let us note that, while v at tree level has the usual relation to the mass 
of the W boson niw, 

mw sin 6w 
V = 2 , 

e 

with 6w being the weak mixing angle and e the electromagnetic coupling, generally does 
not coincide with any of the masses of the physical scalar modes. 

On closer examination of the scalar spectrum, two very different cases can be distinguished 
and are shown in Fig. [T] where |^q| (which controls the couplings to the SM particles) is 
displayed as a function of the mass of the modes^. If rrih lies below rrib (left hand side of 



^ These figures are closely related to the Kallen-Lehmann spectral density of h which is given by p{s) = 

EZoieoYsis-xi). 
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FIG. 1: Compact HEIDI spectra for the two different scenarios referred to in tlie text below. The 
respective values of ruh and mb are marked with vertical lines. 



Fig. H]), we find that the lightest mode, which always lies below rrib as argued above, has 
|(^ol close to one. The other modes form a tower starting at rub with the couplings of this 
"would-be continuum" being smaller by about an order of magnitude and dropping off for 
higher masses. Increasing g"^ pushes the mass of the lowest mode towards zero and enhances 
the couplings of the would-be continuum, while the mass of the lowest mode approaches m/j 
in the limit with \^q \ — )■ 1, and the other modes decouple from the SM. 

If ruh lies above nib (right hand side of Fig. [1]), the situation looks similar for large values 
of g"^. However, if we decrease this parameter, the lightest mode moves towards nib with 
its coupling to the SM decreasing. At the same time, a peak in the would-be continuum 
is appearing at m^. In the limit g"^ — )■ 0, this peak evolves into a single mode with a SM 
coupling of one, while the other modes (this time including the lowest one) decouple. 

For both cases, the asymptotic spacing between the modes is given by the compactification 
scale R~^. In the limit of large R, the modes in the would-be continuum approach each other, 
while their couplings to the SM go to zero. 



C. Considerations regarding simulations 

A thorough study of the phenomenology of compact HEIDI requires its implementation 
into an event generator and the simulation of the different Higgs production channels within 
the HEIDI scenario. Considering that all relevant processes feature at least a four-particle 
final state already at the parton level and taking into account that the number of Feynman 
diagrams is multiplied by the number of propagating scalar fields, speed is a critical issue 
here. Therefore, as the parton-level event generator Whizard 0, 0] can handle final states 
with six to eight particles in the SM and in the Minimal Supersymmetric Standard Model 
(MSSM) [8-10|, and has been used especially for many alternative models of EWSB 11-14 



it is a good choice for the investigation of compact HEIDI models as well. 

For the actual implementation of the model, the tower of scalar fields must be truncated 
at a cutoff scale A. However, if we choose the bulk mass nib to be of the order of the 
weak scale and is small, then it is easy to see from Eq. f lTOj) that the spectrum will 
contain a large number of scalar fields below A. As a consequence, even though at the 
Lagrangian level the HEIDI model is a straightforward deformation of the scalar sector of 
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the SM, consisting of the simple addition of a single five-dimensional scalar singlet, the 
implementation of the model into Whizard (as well as into any other matrix element 
generator) will result in the addition of a very large number of particles and interactions to 
the existing SM implementation of the code. In addition, the number of fields below the 
cutoff depends on the parameters of the model. We therefore choose to follow the approach 
introduced in Ref. [l6| and implement the model using the newly developed FeynRules 



interface to Whizard 17j. In the next section we briefly outline some of the details of the 



implementation using that interface. 



III. PREREQUESITES FOR SIMULATION 

In this section we describe the implementation of the compact HEIDI model of Section 
iniinto FeynRules. We first explain how the numeric spectrum and widths of the scalar 
fields are obtained for a given point in parameter space, while the actual FeynRules imple- 
mentation is presented afterwards. The FeynRules model files are available for download 
from the FeynRules model database located at 

https : //server06 . f ynu.ucl . ac .be/pro jects/feynrules/wiki/ModelFiles 



A. Spectrum and widths 

According to Eq. IQ, the masses of the physical scalar modes are encoded as the 
zeroes of a transcendental function /(A^). As no analytic expression for these roots exists, 
a procedure for their numerical determination is necessary. 

As has already been observed in Section III Al the lowest mode is always located in the 
interval between and rub, and the would-be continuum which starts at nib has exactly one 
mode in each of the intervals of Eq. ( !T0|) . As /(A^) is monotonous between the boundaries of 
the intervals, it is possible to calculate the masses by the bisection method: for each mode, 
we start in the middle of the respective interval and then, depending on the sign of /, split 
the interval we know to contain a root in two until sufficient precision has been reached. 
This is guaranteed to find all modes, starting from the lowest one. 

Once the mass of a mode A^ has been determined, the mixing parameters defined in 
Eq. (ITT]) can be calculated as functions of A^ through application of Eq. (lAlip and Eq. f lAlSp . 
From the mixing parameters, the trilinear and quartic scalar self-couplings can be derived 
according to the Feynman rules given in Appendix IA3[ We have implemented the calculation 
of the masses, mixings and self-couplings into a spectrum calculator. While this is included in 
the model file for FeynRules, it is also available as a stand-alone Mathematica package. 

In addition to masses and couplings, the widths of the scalar modes are required for sim- 
ulating collider cross sections. As the couplings of the (/)k to the SM particles are identical to 
the respective Higgs couplings up to factors of (c/. Appendix IA 3p . their decay widths into 
SM particles can be approximated^ by the corresponding Higgs widths, scaled by factors of 
1^0 I ^- We have chosen to take the necessary SM Higgs width from Hdecay fisl]. In order 
to avoid a run-time dependence on the Hdecay package, we have parametrized the Higgs 



^ There is a small correction coming from diagrams involving both SM particles and other scalar fields ipk- 
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mass dependence of its output in a dedicated Mathematica module (which is also con- 
tained within the FeynRules model file). As the presence of thresholds spoils any simple 
ansatz for fitting, we parameterize the width by a grid of 100 points for Higgs masses below 
500 GeV (with linear interpolation between the grid points), while an interpolating polyno- 
mial of degree 6 is used above 500 GeV. Comparison shows that the result approximates the 
corresponding Hdecay output well at least up to 4 TeV. 

Using this parameterization, our spectrum calculator is able to automatically determine 
the widths of the HEIDI scalar fields. In addition to the aforementioned decays into SM 
particles, two body decays into two scalar modes are also included into the calculation. 



B. FeynRules implementation 

Since the multi-scalar couplings are only of small phenomenological relevance and at the 
same time slow down both FeynRules and the matrix element generator considerably (their 
number growing combinatorically with the number of modes), we restrict the implementation 
to the SM sector of the model. To this end, we insert Eq. f lT2|) into the SM Lagrangian and 
induce in this way the coupling between the HEIDI modes (pk and the SM fermions and 
vector bosons, which is enough to study most of the LHC phenomenology of this model. In 
addition, we restrict the implementation to unitarity gauge. In the rest of this section we 
describe the FeynRules implementation of the model which is special in the sense that the 
number of modes below the cutoff and thus the number of particles in the model is variable. 

Since the compact HEIDI model consists of the SM to which we add the HEIDI sector, 
the FeynRules implementation can be achieved in a natural way by extending the already 
existing implementation of the SM in FeynRules, contained in the file SM.fr included in 
the FeynRules distribution. In this process^, SM.fr can remain unchanged, and the new 
model file only has to encode the new parameters, fields and the modified SM Lagrangian. 

Before the spectrum of the model can be calculated, the Higgs vacuum expectation value 
V and its pre- mixing mass rrih, the bulk mass ruf,, the HEIDI coupling g'^ and the com- 
pactification scale must be specified (c/. Section III Ap . In addition, for the actual 
implementation, we need to truncate the infinite tower of Kaluza-Klein states to a finite 
number by choosing a cutoff scale A. These six parameters, which must be known before 
any code for a matrix element generator can be emitted, have been implemented as a set 
of flags which must be specified before the model can be loaded successfully and which are 
listed in Tab. HI It is important to realize that those values are hardcoded in the generated 
model files. In particular, the electroweak scale implied by the runtime input parameters 
must match the value of v set by the aforementioned flags. 

After setting the flags from Tab. [H the model such defined can be loaded into FeynRules 



and the Whizard interface can be invoked as described in 17 . 



^ As the SM Higgs field h is no longer a mass eigenstate in the model, it can be removed from the set of 
physical fields by adding Unphysical -> True to its definition in SM.fr. However, this step is strictly 
optional and no harm is done by omitting it. 
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HEIDI$cs 


compactification scale R^^ 


HEIDI$mb 


bulk mass 


HEIDI$g2 


5D trilinear coupling squared 


HEIDI$cutoff 


cutoff scale A 


HEIDI$nmodes 


Number of HEIDI modes below the cutoff scale A (determined 
automatically if HEIDI$cutoff is set) 



TABLE I: The flags in the FeynRules implementation of the compact HEIDI model which must 
be set before the model can be used. 
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FIG. 2: Diagrams contributing to signal and background for scalar particle production in pp — ) 
ZZ — )• 41. The shaded box represents the effective gluon fusion operator obtained from integrating 
out the top quark. 



IV. HEIDI SCALARS IN THE GOLDEN CHANNEL AT THE LHC 

In order to study the impact of HEIDI on one of the major Higgs search channels at the 
LHC, we have performed simulations of scalar particle production in gluon fusion with the 
subsequent decay into four leptons via two virtual Z bosons, the so-called golden channel 
for Higgs discovery at the LHC. As gluon fusion is a loop induced process with the dom- 
inant contribution at leading order coming from the top quark triangle, we have added to 
the Lagrangian the corresponding operator obtained by integrating out the top quark (c/. 
Appendix El]). 

At lowest order in the strong coupling, the signal in this process is mediated by gluon 
initiated s-channel type diagrams such as the one shown in Fig. [2] left, while the background 
processes always have quarks in the initial state and consist of diagrams such as those in 
Fig. 12] right. Therefore, a significant part of the background is contributed by processes with 
a valence and a sea quark in the initial state. As those processes tend to be more strongly 
boosted when compared to processes initiated by sea quarks or by gluons, it is possible to 
remove part of the background by a cut on the total rapidity. We thus chose to apply a cut 
of 

bl <2 
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FIG. 3: The scalar field spectra at the four benchmark points at which simulations have been 
performed. 



on the total rapidity of the final state. 

Additional background can be removed by cutting the invariant mass of the leptons to 
the Z mass window. However, as we do not want to lose sensitivity to the scalar modes 
which lie below the pair production threshold, we must allow for lepton pairs with invariant 
mass significantly below the Z mass and therefore have enforced the cut 

10 GeV < mil < 100 GeV 

on the invariant mass of the lepton pairs. If the final state consists of four leptons of the same 
generation, we demanded that at least one of the two possible combinations satisfies the cut. 
In order to account for detector acceptance and for the fact that we require four separately 
resolved leptons, we have applied a px cut as well as cuts on the angular separation and on 
the polar angle of the leptons 



Pt > 5 GeV 



I cos I < 0.99 



Icos^/I < 0.99 



We have performed simulations at four different benchmark points in parameter space 
which are shown in Fig. |3l The parameter sets la and lb differ only by the compactification 
scale and have nih = 200 GeV and = 120 GeV. As discussed in Section fllBI this implies 
that the would-be SM Higgs boson corresponds to a bump in the would-be continuum around 
ruh as clearly visible in the spectra shown in Fig. 12] left. At point la, the compactification 
scale is = 10 GeV with more than 20 modes lying below 250 GeV, while we only have 
6 such modes at lb with = 50 GeV, their couplings to the SM particles being enhanced 
in exchange. 

At Ila and lib, we set rrih = rub = 160 GeV. As discussed in Section HTB] and clearly visible 
from Fig. [3] on the right, the spectrum is very different in this case, with the lowest mode 
corresponding to the would-be Higgs boson, and the couplings of the would-be continuum to 
the SM being strongly suppressed. Again, we chose = 10 GeV for Ha and = 50 GeV 
for lib. 

For the actual simulation, we took the first two quark generations and gluons into account 
as potential initial particles, simulated the full 2 — )■ 4 processes using the CTEQ6L [l^ 
series of parton distribution functions and generated unweighted events for an integrated 
luminosity of 100 fb~^ at ^/s = 14 TeV. As the background is the same for all four benchmark 
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FIG. 4: Total invariant mass distributions in pp — )• 41 at benchmark points la and lb. For the 
distributions in the right column, a Gaussian smearing with a = 2 GeV has been applied to the 
invariant mass in order to simulate measurement errors. 



points, we generated a single set of background events, while a fresh set of signal events was 
generated for each point. 

The top left panel of Fig. H] shows the resulting invariant mass distribution for benchmark 
point la. Despite the fact that the scalar field tower already starts at 120 GeV, the actually 
visible resonances only start at the Z pair production threshold, demonstrating that the 
SM couplings of those modes are much too small to win over the suppression coming from 
the off-shell Z propagators. Above the threshold, the comblike structure coming from the 
scalar resonances is clearly visible not only in the actual signal, but persists also when the 
background is added. 

The same distributions are shown in the top right panel of Fig. HJ but now with a Gaussian 
smearing with standard deviation a = 2 GeV applied to the invariant mass in order to 
simulate a measurement error. Although the width of the Gaussians is still significantly 
smaller than the average mode spacing of about 10 GeV, the histogram shows that this is 
already enough to destroy the comblike structure of the would-be continuum. What remains 
is a broad excess above the background which peaks around nih and resembles a single, very 
broad resonance. 

The bottom panels of Fig. Hlshow the unsmeared (left) and smeared (right) distributions 
for benchmark point lb. As detailed above, the only difference with point la is the larger 
spacing of the individual modes and the therefore enhanced couplings to the SM. Before 
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FIG. 5: Like Fig. [U but for the benchmark points Ila and lib. 



the smearing is applied, four modes are visible in the signal (one of them below the Z pair 
threshold with a mass of ~ 150 GeV), and the three heaviest of them also stick out above 
the background. While the smearing distorts the shape of the peaks and makes them less 
pronounced, the smeared peaks are still visible above the background. 

In the same way, the invariant mass distributions for benchmark points Ila and lib are 
shown in Fig. [5l the top panels again showing the case with the smaller mode spacing. The 
different structure of the spectrum discussed above also manifests itself in the distributions. 
The lowest mode, which again lies at about 120 GeV, has nearly SM-like couplings in this 
scenario and is therefore a prominent feature in the signal distribution. Despite lying far 
below the Z pair threshold, the corresponding peak is more than comparable in size with 
the peaks from the would-be continuum. After adding the background, the lowest mode still 
sticks out and remains visible, while the comblike structure of the other modes is masked 
by the background fluctuations with only a diffuse excess being observable. Applying the 
smearing degenerates the peak of the lowest mode and homogenizes the (small) excess over 
the background. 

The distributions for benchmark point lib again differ by the mode count. While the 
appearance of the lowest mode is virtually identical to scenario Ila (as is clear from Fig. |3]), 
the would-be continuum now exhibits two moderate peaks in the signal, both of which 
survive when the background is added. After the smearing is added, at least one of those 
two peaks remains visible over the background together with the lowest mode. 
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Our simulations show that the model could, depending on the scenario realized, appear 
in very different guises, ranging from a Higgs-like peak with reduced magnitude and a diffuse 
continuum (Ila), a series of sharp resonances (lib and lb) or a single, very broad excess (la). 
At least for la / lb, the signal is clearly visible over the background and a discovery might 
be possible in this channel. For lib and especially Ila, the situation looks more pessimistic; 
especially in Ila where the diffuse shape of the signal might get lost in the background. 
However, while being the first detailed exploration of a HEIDI-like scalar sector at the LHC, 
our simulations are still too simplistic to make realistic efficiency estimates for a discovery. 
To this end, a more detailed study of the background and of higher order corrections is 
required as well as the inclusion of other potential discovery channels like gauge boson 
fusion and "HEIDI strahlung" . We postpone this to a future publication. 

V. CONCLUSIONS 

In this paper, we have presented a first study of the LHC phenomenology of an extended 
scalar sector built from a compact extra dimension, where the SM fields are living on a 
brane and only a single singlet scalar field is allowed to propagate in the bulk. This is a 
representative of a class of models known as higher-dimensional ( "HEIDI" ) Higgs models. 

Despite the presence of just a single higher-dimensional scalar singlet, the physics and 
phenomenology of this model is astonishingly rich. The main effect comes from the Kaluza- 
Klein tower of scalar resonances mixing with the SM Higgs scalar field. As the extra- 
dimensional scalar fields couple to SM particles through this mixing, the SM Higgs boson 
gets diluted to a tower of scalar fields of ascending mass whose couplings to the SM particles 
are reduced significantly compared to the SM Higgs boson. 

We concentrated on the so-called "golden channel" for a not too light Higgs boson decay- 
ing predominantly into two Z bosons, the latter decaying into four leptons, and restricted 
ourselves to the gluon-fusion production channel. Our main results in Section |IV] show a 
variety of different shapes which such a model could produce, depending on the point in 
parameter space under consideration. These shapes range from a simple reduction of a SM 
Higgs signal with a diffuse continuum tail to a series of Kaluza-Klein peaks to a very broad 
continuous bump, showing the phenomenological richness of these otherwise simple exten- 
sions of the SM. While this phenomenological variety is a generic property of the model 
which derives from its spectrum, it would certainly be interesting to see how the other clas- 
sic Higgs production channels perform for this type of model. We conclude by remarking 
that even a rather trivial deformation of the SM in the scalar sector can have the potential 
to jeopardize the common assumption that the LHC cannot fail in its mission to discover 
either the Higgs boson or new physics. 
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Appendix A: Diagonalizing compact HEIDI 
1. Vacuum expectation values 

From Eq. ([1]) and Eq. (|3]), we can read off the classical potential for the vacuum expec- 
tation values of the fields Eq. @, 

V{v, w) = ]xv' - lf,'v' + V (Imlwl - -^v'w,) . (Al) 



For any fixed value of v, Eq. f lAip decomposes into a constant term plus a sum of second 
degree polynomials V^, each depending only on a single Wk- If we have rUb > 0, then each of 
the Vk is bounded from below, and the the minima are given by 

Wk = — . (A2) 



Therefore, the potential V is always bounded from below as a function of the Wk for fixed v 
and ruh > 0, and inserting Eq. ( ]A2l) into Eq. ( lAlj) gives the minimum 

Voiv) = \iX-a)v'-l^^'v\ (A3) 

With the help of 

°° 1 1 TTX 

fc=0 ^ + ^ ^ 

a in Eq. (IA3P can be written as 



Q! = — > — ^ — ^ = — coth (Rnmu) . (A5) 



In order to derive the conditions under which Eq. ( lAip is bounded from below as a 
function of both v and the Wk, it is sufficient to examine the asymptotic behavior of 
Eq. (lASp which is determined by the sign of the quartic term. As a can be easily seen to 
be positive definite from Eq. f lASp . we end up with four different scenarios: 



1. > , A > a: 

V is bounded from below and has two minima at 



The point v = = Q is unstable. 
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2. ^2 > , A < a: 

The potential is not bounded. 

3. ^2 < ^ X>a: 

V is bounded from below and has a single minimum at v = Wk = 

4. < , A < a: 

The potential is not bounded. 

Therefore, we need //^ > and A > a in order to have a stable theory with spontanous 
symmetry breaking. The two solutions Eq. flA6l) are equivalent and only differ by a field 
redefinition h — )■ —h as can be seen from Eq. ([7]), and we therefore are free to make the 
choice Eq. ([5]). 

2. Mass eigenstates and wave functions 

Arranging the fields h and Uk into the vector 

T = {h,uo,ui, . . . 
and introducing the mass matrix M 




for k = 
for > 



(AT) 



with the abbreviation 



2gv 



NkVR 



the mass and mixing terms in Eq 



([7]) can be written as 




In order to find the eigenvalues A^ to the eigenvectors 
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(A8) 



we rewrite the eigenvalue condition of Eq. flA8|) in terms of the components ^ 



oo 




(A9) 



j=0 



A^ef = P.-i^o' + for 



i > . 



(AlO) 



^ Note that this definition of the is equivalent to that given in Eq. ([IT]). 
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Since assuming = and ^ for some i > leads to a contradiction, we can therefore 
solve Eq. f lA10|) in terms of 

= ^oTT^ ■ (All) 

Upon inserting Eq. ( lAllI) into Eq. (]A9I) . .^q drops out and application of Eq. ( ]A4[) finally 
leads to the result Eq. ([9]). In order to obtain the normalized eigenvectors (^^, we rewrite the 
normalization condition using Eq. flAlip 

The sums appearing in Eq. (]A12p can be performed analytically, and solving this equation 
for ^Q, one obtains. 

The other components of the eigenvectors can be derived by inserting Eq. (1A13P into 
Eq. (lATT]) . 



3. Feynman rules 

As the model differs from the SM in the scalar sector only, all couplings involving fermions 
and gauge bosons are unchanged and need not be repeated here. To obtain the Feynman 
rules involving the scalar fields, we have to take the full Lagrangian of the model (which 
consists of the SM Lagrangian together with the scalar sector Eq. ([7])) and express h and 
the Uk through the mass eigenstates 0^ by application of Eq. (fT2|) . 

For the couplings between the (pk and the SM vectors and fermions, the resulting 
Feynman rules are trivially obtained by taking the SM rules and replacing the Higgs legs 
with the (j)k, multiplying with a factor for every scalar leg. The results are Feynman 
rules which look like 



SO ■ ■ ■ SO ' 



where the double line represents all vector and fermion lines, and F is the SM vertex factor 
(including the color and Lorentz structures). 

The trilinear scalar couplings are more complicated: the piece of Eq. ([7]) which encodes 
them is 

oo 

3 1,2 




To obtain the corresponding Feynman rule in a compact form, we insert the physical fields 
via Eq. ( 1T2|) . sum the series using Eq. (lAlip and Eq. ( ]A4p . apply Eq. ([9]) and replace A 
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using Eq. ([6]). We eventually derive 



ko^eo (A^ + A? + - 6ml) 



The quartic couplings are considerably simpler to obtain since it is sufficient to insert 
Eq. f[T^ into the /i^ term of Eq. ([7]) . We directly get 




-^3!A aaeoeo ■ 



4. The effective scalar-gluon-gluon coupHng 

Straightforward evaluation of the matrix element for on-shell gluon fusion of a single 
scalar mode (pi via a top quark loop gives (c/. |22|) 

a,ki,iJL a, fci,/x 




b,k2,v b,k2,v 



^il n (1 + (1 ; ^0 9{r.^\ jk'ik^ - 9'" {hk^)) Sg, (A14) 
pin) 

where the function g{T) is defined piecewise as 




for r > 1 

for r < 1 



and Ti is defined as 
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with the top quark mass nit and the scalar field mass Aj. The amplitude Eq. f lA14p is 
reproduced by the trilinear part of the gauge-invariant effective operator 



^99<t> = E T^^oP(^.) <P^ tr G^.G^^ (A15) 
^-^ Ibnv 

i=0 

with the gluon field strength tensor C^^. 

The scalar (pi can be taken off the mass shell simply by replacing by its momentum 
p^, and the lowest order of the expansion of Eq. (lA14p in can then be obtained by taking 
the limit 

2 

lim p{t) = - 

Inserting this limit into Eq. (]A15|) . we obtain the lowest dimension operator contributing to 
gluon fusion in the effective theory obtained by consistently integrating out the top quark 

oo 

^.V = E^?o^^trG,.G- (A16) 

i=0 

Both effective operators are available in our FeynRules HEIDI implementation via the 
function LHEIDIgg. Calling this function as LHEIDIgg ["heavytop"] generates Eq. flA16l) . 
while Eq. (lAlSP can be obtained by omitting the argument. For the simulation results 
presented in this work, Eq. (I A 161) has been used. 
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